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ABSTRACT 


Many satellites, such as ATS-A, ATS-D, DODGE, RAE, etc., have 
clamped to the surface of their central body long elastic members which 
when set into vibratory motion can significantly influence the motion of 
the central body. If such a satellite is reasonably symmetrical, we can 
predict its dynamic characteristics by modeling the entire satellite as 
a symmetric double-beam system; that is, by a rigid symmetric central 
body having clamped to its surface two long, diametrically opposed uni- 
form elastic beams with tip weights possessing identical physical and 
geometrical properties. In this paper, dimensionless equations which 
define the natural modes and frequencies of such a satellite system are 
derived and solved. The solutions are outlined in graphical form and 
then used to solve the equations that describe the elastic response of the 
satellite to an arbitrary periodic forcing function. The results in both 
graphical and analytical form make it possible to predict with slide-rule 
accuracy the natural frequencies of any satellite that can be modeled as 
a symmetric double-beam system. The dynamic response equations 
have been applied to a particular problem associated with ATS-D. It is 
shown that if the ATS-D microthruster were operated at any frequency, 
including the system resonants, the viscoelastic damping of the elastic 
beams would be sufficient to prevent the resonant-oscillation amplitudes 
from exceeding the elastic limits of the attached appendages. 
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THE DYNAMIC CHARACTERISTICS OF SATELLITES 
HAVING LONG ELASTIC MEMBERS 


by 

Harold P. Frisch 

Goddard Space Flight Center 


INTRODUCTION 

The design and computer simulation of an active control system for a satellite having long 
elastic members (booms) requires knowledge of each elastic mode of vibration that can be excited, 
and its associated natural frequency. This is essential for system design, and helps to develop 
and verify the computer simulation of the entire coupled dynamic system. 

To predict with accuracy each vibration mode and resonant frequency, we must account for 
the geometric configuration of the satellite system; the mass, size, and rotational inertia of the 
rigid satellite body; and the elastic properties of the appendages. Frequently the geometrical con- 
figuration permits describing the satellite by a symmetric double-beam model, that is: two dia- 
metrically opposed uniform elastic beams with tip weights, attached to a rigid central body having 
mass, finite size, and rotational inertia (Fig- 
ure 1). It is this type of satellite system that 
will be extensively discussed in this paper. 

The equations developed in dimension- 
less form and the results given provide suf- 
ficient information to permit computing, with 
slide- rule accuracy, the natural frequencies of 
each of the pertinent modes of vibration for any 
such satellite. By comparing the natural fre- 
quencies with those resulting from the limiting cases of zero or infinite satellite body weight, ro- 
tational inertia, or tip weight, it is possible to justify further simplifications in a detailed com- 
puter simulation of the entire satellite control system. 

Such a system (as shown in Figure 1) is symmetrical with respect to a plane passing through 
its center of mass. Therefore, the natural modes of vibration can be separated into two groups: 
those that are symmetric about the axis of symmetry, and those that are antisymmetric about the 
axis of symmetry relative to the point on the axis that is coincident with the system's center of 
mass. The axis of symmetry is defined as that which is normal to the undeflected shape of the 
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Figure 1— Symmetric double-beam satellite system. 
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system, passing through the center of mass and in the plane of vibration. This symmetry permits 
one to derive the two frequency equations and the mode shapes of the entire system by simply solv- 
ing the vibration equation over the elastic portion of one beam twice: once with the boundary con- 
ditions associated with the symmetric modes of vibration, and again with the boundary conditions 
associated with the antisymmetric modes of vibration. 

The derived modes form a set of orthogonal functions; these can be used to obtain a series 
solution to the equation describing the vibrational motion of each point of the system. 

To obtain the solution the applied force distribution is decomposed into an equivalent force 
system. That is: 

1. A force distribution, symmetric about the system T s axis of symmetry, which excites only 
the symmetric modes of vibration, and 

2. A force distribution, antisymmetric about the system T s axis of symmetry, which excites 
only the antisymmetric modes of vibration. 

It follows from the principle of superposition that the system’s dynamic response to the equiv- 
alent force system is identical with its dynamic response to the actual force distribution. 

The derived equations, to illustrate these results, are used to solve a particular problem as- 
sociated with ATS-D. They are specialized so as to describe the response of the satellite system 
to a microthruster of periodic intensity located at the point where one of the booms is attached to 
the rigid central body. It is shown that the viscoelastic damping of the boom material itself is 
enough to prevent the peak bending moment at every point of the system from exceeding the critical 
buckling moment of the boom itself, even if the thruster excites a system resonant frequency. 

1. BASIC THEORY 

The equations and results presented here can be applied to any satellite that can be modeled 
as the symmetric double- beam system shown in Figure 1. This is an elastic system, which has 
two identical diametrically opposed uniform elastic beams with tip weights, clamped to the surface 
of a rigid central body. The central body is assumed to be perfectly rigid and have finite size, 
mass, and rotational inertia. Thus the elastic system under study is treated as one having dis- 
continuities in mass and stiffness along its length. 

The present analysis first derives the orthogonal functions that define the normal modes of 
vibration and the corresponding frequency equations, in a form that may be readily applied to any 
similar system. These equations are then used to analytically solve the damped forced- vibration 
equation. The partial differential equation for the damped vibrational motion of every point along 
the systems length, when it is excited by a known force distribution, is 


d 2 <9 2 Y(x, t ) (9 2 Y(x , t ) dY(x, t) 

2 EI ( X ) ZT2 + ^( X > TT2 + £( x ) a? ” F ( x ’ * > 


<9x 


dx 2 


<9t 2 


( 1 . 1 ) 
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where 


X 

t 

Y(x, t) 
F(x, t) 
P(x) 
EI(x) 
/5(x) 


position coordinate measured positive from the boom tip, ft (Figure 1) 
time, sec 

transverse displacement at (x, t), ft 
external force distribution at (x, t), lb ft" 1 
mass distribution at x, slugs ft -1 
bending stiffness at x, lb ft 2 
viscous damping coefficient at x, lb sec ft" 2 .* 


When the viscous damping coefficient /3(x) is proportional to the mass distribution p(x) equation 
1.1 can be solved by separation of variables; the resulting solution is of the form 


Y ( x, t) = 2> (t)y n ( >0 • ( J -2) 

n = 1 


where 

y n (x) is the nth normalized mode of the undamped unforced equations of motion of the entire 
system, and 

a n (t) is the time-dependent generalized displacement coordinate associated with the nth mode, 
feet. 

These mode shapes y n (x) that are to be used in Equation 1.2 must be derived from the equation 


j* EI(X) 


d 2 Y(x, t) 
ax 2 


* «*> 


(1.3) 


and must satisfy the time- independent boundary conditions of the system. 

Because, when the system vibrates in any one of its normal modes of vibration the deflection 
of every point of the system varies harmonically in time, the solution to Equation 1.3 is of the 
form 


Y(x, t ) - y(x) e* wt , 


(1.4) 


where y(x) defines the shape of the normal mode and 


co is its frequency of vibration (radians sec' 1 ). 


*For a complete symbol list, see Appendix A. 
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Directly substituting Equation 1.4 in Equation 1.3 gives the differential equation 


EI(X) 


d 2 y(x) 
dx 2 


- p(x) <x> 2 y(x) 


(1.5) 


Equation 1.5, along with its time- independent boundary conditions, defines an eigenvalue prob- 
lem. The eigenvalues, which are the solutions to the characteristic (frequency) equation of the 
problem, define the natural frequencies of the system, and the eigenvectors associated with each 
of the eigenvalues define the orthogonal modes of the system. 

Since Equation 1.5 has parametric discontinuities along its length which cannot be assumed 
negligible a geometric construction must be used to develop solutions to this equation. 

Reference 1 states that the mode shapes that are assumed to satisfy the equation 


d 2 d 2 y n (x) 

d ^ 2 EI(x) ~~dT 2 " y n ( x > = 0 


( 1 . 6 ) 


must satisfy two basic criteria: they must impose as little constraint as possible on the motion 
of the system, and they must satisfy the kinematic boundary conditions. It is also extremely use- 
ful but not essential that the constructed modes be orthogonal. 

The actual mode shapes y n (x), defined over the entire system length, that satisfy Equation 1.6 
are obtained as follows: 

1. Note that the complete set of natural modes of vibration can be divided into two groups: 

a. Those symmetric about the system’s axis of symmetry: "symmetric modes of 
vibration." 

b. Those antisymmetric about the system’s axis of symmetry relative to the point coinci- 
dent with the system’s center of mass: "antisymmetric modes of vibration." 

2. Let 

El = bending stiffness, uniform along length of elastic beam, lb ft 2 , 

P = mass distribution, uniform along length of elastic beam, slugs ft** 1 , 

then solve the eigenvalue problem defined by the equation 


El 


d 4 y(x) 
dx 4 


pa) 2 y(x) 


0 


(1.7) 


4 


over the elastic portion of half of the system; that is, over the length of one beam twice, as 
follows: 

a. Use the boundary conditions that define the effect of the central body T s mass on the 
beam at its point of attachment, and the effect of the tip weight on the beam at its 
point of attachment, when the system is vibrating in one of the symmetric modes of 
vibration. 

b. Use the boundary conditions that define the effect of the central body’s size and rota- 
tional inertia on the beam at its point of attachment, and the effect of the tip weight on 
the beam at its point of attachment, when the system is vibrating in one of the anti- 
symmetric modes of vibration. 

3. Reflect about the system 1 s axis of symmetry each derived eigenvector of Equation 1.7, 
symmetrically or antisymmetrically as dictated by the boundary conditions. The two points anal- 
ogous to the points at which the beams are attached to the central body will then be displaced by 
a distance 2R, the length of the central body, from each other. Connecting these two points with a 
linear segment generates a function that is definable at every point along the system’s length. 

This function satisfies Equation 1.6 at every point along the flexible portion of the system; it 
satisfies all boundary conditions and accounts for the length and rigidity of the satellite body. 
Furthermore, since the eigenvectors of Equation 1.7 are orthogonal and the satellite length is 
small compared to the system length, the symmetric and antisymmetric modes will form two sets 
of orthogonal functions. It is easily shown that these two sets are mutually orthogonal by noting 
that the integral of the product of an odd and an even function over symmetric limits is zero. 

Note that each mode shape as constructed is not unique but defined only up to a multiplicative 
constant. The mode shapes, however, can be made unique by requiring that the orthonormality 
condition 



/>(*) y n O) y m O) dx 


m t s 

T n,m 


( 1 . 8 ) 


be satisfied for all integers n and m, where 


m t = total mass of system, slugs, 

2L = total length of system, ft, 

ro if n / m, the Kronecker delta function 
n,m 1l if n = m. 


5 



2. NATURAL MODES AND FREQUENCIES 


As previously mentioned, the natural modes and frequencies of the system are constructed 
from the solutions to Equation 1.7 along with the appropriate boundary conditions that define the 
symmetric and antisymmetric modes. 
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In order to avoid notational difficulties it 
is convenient to define the parameter z over the 
elastic portion of half the system, as shown in Fig- 
ure 2. Since z is linearly related to x, Equa- 
tion 1.7 can be rewritten as 

£i d 4 y*(z) _ pa)2yNc(z) = Q (2.1) 

dz 


Figure 2— Coordinate system used to derive 
frequency equations. 


and the boundary conditions defined at the points 
z = 0 and z = L 0 = L - R , L 0 being the length of 
one beam. This equation has the general solution 


y* (z) - Cj cosh ^z + C 2 sinh kz + C 3 cos kz + C 4 sin \z , 


( 2 . 2 ) 


where 


\ 4 


poj 2 

la" 


(2.3) 


and C l9 C 2 , C 3 , C 4 are constants of integration. 

The actual values of k to be used in Equations 2.2 and 2.3 are the solutions to the frequency 
equations of the system, which are derivable from the equations defining the boundary conditions. 

The boundary conditions for the two cases mentioned are: 

Symmetric modes of vibration: 

dy* ( z ) _ n 

dz " 0 ’ 

z = 0 


- El 


d 3 y* (z) 
dz 3 


W s 

+ “5g y* ( z ) 


o , 


EI 


d 2 y* ( z) 
dz 2 


= 0 , 


EI 


d 3 y* (z) W T 


+ — y* (z)J 

z=L n 


= 0 


(2.4) 


6 



Antisymmetric modes of vibration: 


dy* (z) 
dz 


1 

r y*( z ) 


z = 0 


o , 


- El R 


1 y* (z) 

dz 3 


+ El 


d 2 y* ( z) 
dz 2 


I 

s 

2R 


P* ( z )| 

z = 0 


El 


d 2 y* (z) 


dz 2 


0 , 


El 


d 3 y* ( z ) 
dz 3 



oj 2 y* ( z ) 


0 , 


(2.5) 


where 


L 


o 


W T = tip weight, lb, 

W s = weight of rigid satellite body (without beams and tip weights), lb, 

l s = mass moment of inertia of rigid satellite body (without beams and tip weights) 
about the axis through the center of mass and normal to the plane of vibration, 
slugs ft 2 , 

g = acceleration of gravity, ft sec' 2 , 

L - R = length of one beam, ft, 

2R = length of satellite body, ft, 

2L = length of entire system, ft. 


The frequency equation associated with each case is derived by substituting the general solu- 
tion (Equation 2.2) in the four equations defining the boundary conditions. This yields four homo- 
geneous equations in the four unknowns, C lf C 2 , C 3 , C 4 , for each case. It follows that for a non- 
trivial solution to exist, the determinant of the coefficients of C 1? C 2 , C 3 , C 4 must be identically 
equal to zero. The evaluation of this determinant yields a transcendental equation in commonly 
referred to as the "frequency equation." Thus: 


Frequency equation for symmetric modes of vibration: 

A direct substitution of the general solution, Equation 2.2, in the four boundary conditions, 
Equations 2.4, yields: 
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0 


c + c = 

C x coshXL 0 + C 2 sinh\L 0 - C 3 cos ^L 0 ~ C 4 sin\L 0 = 0 

W 

- EIX 3 (C,-C 4 ) + 21^(0, + C 3 ) = 0 

EIX 3 sinh XL 0 + C 2 cosh XL„ + C 3 sin XL 0 - C 4 cos XL 0 ) 

W T 

4- ~ ~ co 2 (Cj cosh XL 0 + C 2 sinh XL 0 + C 3 cos XL 0 + C 4 sin XL„) = 0. (2.6) 

For a non- trivial solution to exist, the determinant of the coefficients must be equal to zero; 
that is 


<lii 
^21 
<*31 
^4 1 

where 

= 0 

q i2 = 1 

<ll3 = 0 

^14 = 1 

^31 ■ 2W b XL o 


^33 “ 2Wg ^ L o 

<*34 = 1 


<*12 

^1 3 

^14 

^22 

^2 3 

q 24 

^3 2 

^33 

^34 

^42 

%3 

q 4 4 


q 2 J - cosh XL 0 

q 22 = sinhXL 0 

q 23 = - cos XL 0 

q 24 = -sin XL 0 

W T 

q 41 - sinhXL 0 + \fj XL 0 cosh XL 0 

i. 

q 42 - cosh XL 0 + XL 0 sinh XL 0 

«r 

q 43 a sin XL 0 + XL 0 cos XL 0 

Wt 

q 44 - - cos XL 0 + w XL 0 sinXL 0 


W fi ' pgL 0 - weight of one beam, lb , 


(2.7) 


X 4 


pod 

TT 


( 2 . 8 ) 
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Directly substituting Equations 2.8 in Equation 2.7 gives the following: 


tanh XL q + tan\L Q + 2W (l + sec ^-L 0 sech A.L Q ) 

W s “ ~~ 

2X.L q - 2 W~ K) 2 (tan \L 0 - tanh \L 0 ) 


This is the frequency equation for the beam of length L 0 having the prescribed boundary conditions. 
Furthermore, because of the symmetry of the system, it is also the frequency equation for the 
symmetric modes of the entire system. 

Frequency equation for antisymmetric modes of vibration: 

A direct substitution of the general solution, Equation 2.2, in the four boundary conditions, 
Equations 2.5, yields: 


M C 2 + C 4 ) " R ( C i + C 3 ) = 0 


C x cosh XL 0 + C 2 sinh\L Q - C 3 cos XL 0 - C 4 sinXL 0 - 0 


-EIR \ 3 (C 2 -C 4 ) + EJA 2 (C, -C 3 ) + -£a> a (C, +C 3 ) = 0 


EIX 3 (Cj sinh \L 0 + C 2 cosh \L 0 + C 3 sin XL 0 - C 4 cos XL 0 ) 


i. ,2 (Cj cosh XL 0 + C 2 s inh XL 0 + C 3 cos XL 0 + C 4 s in \L q ) - 0 


( 2 . 10 ) 


For a non-trivial solution to exist, the determinant of the coefficients must be equal to zero; 
that is: 


( 2 . 11 ) 


I 



where 




R* 


q 21 - cosh \L 0 


2 ~ Q 2 2 " sinh AJL 0 


1 

R* 


q 2 3 


cos A.L n 


1l4 = XL 0 


q 24 - - sin A.L 0 


^31 R* + 2W_ (^o) 2 


q 3 2 ^0 


1 W s 

^33 = "R*" + 2W I (^ L o) 2 


^3 4 ' , + ^ L 0 


q 41 - sinh\L 0 + yj— A.L 0 cosh\L 0 


q 42 - cosh \L Q + ^ ^L 0 sinh\L 0 


q 43 - sin \L 0 + Wb Ai 0 cos\L 0 


q 44 - " cosAL o + Wb \L„ sin \L 0 


W B = PgL 0 3 weight of one beam, lb 


* ^ S S 

W - — T - effective inertia weight of rigid satellite body, lb 

s R 2 


R* - — n normalized satellite radius 


(2.12) 
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Directly substituting Equations 2.12 in Equation 2.11 gives the following: 


Let 


(tanh \L Q - tan \L Q - 2 XL q R* tan XL 0 tanh \L 0 - (XL 0 R*) 2 (tanh \L Q + tan XL Q ) 


2W 1 


y R* 2 (\L 0 ) 3 (l + sec \L 0 sech \L 0 ) 




(2.13) 


2 tan XL 0 tanh XL 0 + 2R* XL q (tanh\L Q + tan XL Q ) + 2(R* \L 0 ) 2 


2W r 


R* 2 (XL Q ) 3 (tanh \L 0 - tan \L 0 ) 


I B = mass moment of inertia of a beam of length 2L 0 about an axis normal to the length and 
passing through its centroid, slugs ft 2 . 


Then 


2 

TPL 0 3 


(2.14) 


This expression can be used to rewrite Equation 2.13 in terms of the inertia ratio l s /l B . That is, 
by substituting 


W * 

s 

“w: 


i s g = K 

pgL 0 R 2 pL 0 3 R* 2 


21 


3I b R* 2 


(2.15) 


(2.16) 


in Equation 2.13: 

(tanh XL 0 - tan \L 0 - 2XL 0 R* tan XL Q tanh XL Q - (\L Q R*) 2 (tanh XL Q + tan XL Q ) 

I s 

" 31^ ( XL o) 3 (l+ secXL 0 sech\L 0 )^ 

B JXL 0 2 tan XL 0 tanh XL 0 + 2XL 0 R* (tanh Xh Q + tan \L 0 ) + 2(XL Q R*) 2 

+ 3 J - (XL 0 ) 3 (tanh\L 0 - tan XL 0 

L B 

This is the frequency equation for the beam of length L 0 having the prescribed boundary conditions. 
Furthermore, because of the symmetry of the system, it is also the frequency equation for the 
antisymmetric modes of the entire system. 


Equations 2.9 and 2.16 are transcendental equations that cannot be readily simplified except 
for extreme values of the dimensionless ratios included. For any particular set of parameters, 
however, the equations can be solved accurately by numerical techniques. 
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It is convenient to put those values of kh 0 that satisfy Equation 2.9 and 2.16, in a particular 
order. That is, let: 

(x_ L o ) 2 n~ i = the nth value of \L 0 (in increasing order of magnitude) that satisfies the frequency 
equation for the symmetric modes of vibration. Equation 2.9, 

(XL 0 ) 2n = the nth value of \l 0 (in increasing order of magnitude) that satisfies the frequency 
equation for the antisymmetric modes of vibration, Equation 2.16. 

Thus, from Equation 2.3, the mth natural frequency of the system is given by 



where it is understood that 


k 


m 


KoL 

Lo 


(2.17) 


(2.18) 


Associated with each eigenvalue k m there exists an eigenvector y m *(z) that satisfies the differ- 
ential equation 


d 4 y * (z) 

-174- ~ K y m * ( z ) = 0 (2.i9) 

and the four boundary conditions given by Equations 2.4 or 2.5. This function y m * (z), along with the 
orthonormality condition, Equation 1.8, is used to construct y m (x), the mth natural mode of the 
entire system. 

The equation defining y m * (z) is derived as follows: 
m odd; symmetric modes: 

Let c lm , C 2m , C 3m , C 4m be the constants of integration associated with the mth natural frequency 
ai m . Directly substituting X m into the first three of Equations 2.6 and solving for the constants leads 
to: 


c = - c 

L 2m 4m ’ 

C lm COshX m L 0 + C 2n, ( sinhX m L 0 +sLn K L o) = C 3m COsX m L 0- 

w s 

2C a m - 2W7 L o ( c i m + C 3m ) = 0 • (2.20) 
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Hence, 


C 2m “ + 4W_ L 0 [ l +G m) C lm - 


C, = G C, , 

3 m m lm ’ 


4W_ 


, L o (1 +G J C lB , 


( 2 . 21 ) 


where 


secX „ L o l 1 + 4 W b ^ m L 0 tanh\.L 


W s 

sechX m L o l 1 " 4Wl \, L o tan ^~ L ‘ 


\ h_ 

oj + 4W B 


X L n sech X L n tan X L n 

m 0 m 0 m 0 


X L n sec X L n tanhX L n 


( 2 . 22 ) 


Substituting these equations into the general solution gives the equation defining y m * (z) as 


y * (z) n C. cosh X z 

m v 7 lm m 


1 + au] X L n (l + G ) t anh X z+G cosX zs ech X z 

t tTVg m u \ m/ m m m m 

w 


4 W ‘X L n (l + G } s in X z sech X z 

t"VY g m U \ m / m m 


(2.23) 


for all odd values of m. The evaluation of the magnitude of C lm will be discussed in Section 3. 
m even; antisymmetric modes: 

Let c im> C 2 m 9 C 3 m> C 4 m constants of integration associated with the mth natural frequency 

oi m . Directly substituting X m into the first three of Equations 2.10 and solving for the constants 
leads to: 


( c i m +c 3m ) + \,l 0 r* (C 2m + C 4 J = 0 , 


w * 


W* 

1 + 2^ R * (^L 0 ) 2 _ 

C , m - \,L 0 R * C 2 m + 

" 1 + 2W^ R * (^m L o) 2 


G 3m + L 0 R * G 4m 


C ln, COsh L 0 + C 2m Sinh \n L 0 " C 3 m COS L 0 ~ C 4 m sin \n L 0 “ 0 ' 


(2.24) 
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Hence, 


where 


" 2 m 1 m 


2 + 2 t R * 2 (^l 0 ) 3 H m 


R *^- L o (2 - ^K*(^ m Lo) 2 ) 


"3m 1 m 


2W, 


R* (k L„) 2 + 2A. L n R* H 

\ m U / m U it 


2 2W_ R * (^m L o) 2 


C A = C. H , 

4m 1 m m 


(2.25) 


R* = 


R* L n 


J o 

K g 

R 2 

f 2secX : L o ( R * \»l 0 +tanh L o) 

W * 

2Wg R * 2 K L o ) 3 (secX m L 0 + sec h ^L 0 )^ 

W* 

2sech\ m L 0 (R^^Lg+tan^L,,) - ^ R* (\ m L 0 ) 2 (tanh \ m L 0 sec \ m L Q 

+ tan\ L n sech A L n ) 

m O m 0 / 


•(2.26) 


Substituting these equations into the general solution gives the equation defining y m * (z) as 


y m *( Z ) = C lm cosh \n z + 


2 + it R * 2 (\» L o) 3H . 


R *\, L o (2 - it R *(\.L 0 ) 


tanh A z 


2 1 R * (\n L o) 2 + 2A m L 0 R*H„ 
“ ' W * 

2 - it R * K L o) 2 


cos A z sech \ z + H sin A. zsechA z> (2.27) 

m m m f 


for all even values of m. 
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3. CONSTRUCTION AND NORMALIZATION OF MODES 


In the preceding section the functions y m * (z), which are to be used to construct the symmetric 
and antisymmetric modes, have been derived up to the multiplicative constant C lm . By following 
the rules outlined for the construction of y m (x), the mth natural mode of the system, and applying 
the orthonormalization condition, Equation 1.8, we can uniquely define the mth mode. 

Figures 1 and 2 and the boundary conditions previously given show that the function y m * (z), 
defined at every point in the interval 


0 < z < L 0 , 

is identical with the function y m (x) at every point in the interval 

L 0 + 2R < X < 2L , 


when 


X - L 0 - 2R . 


The determination of y (x) over the full interval 


0 < x < 2L 

for the symmetric and the antisymmetric modes uses the following construction procedure: 
m odd, symmetric modes: 


For this case, the function y m * (z) is defined 

by Equation 2.23 and 



1 

o 

Vj* 

for 

o 

V| 

X 

V| 

o 

y m ( x > 

- < 

1 

y m * (°) ■ 

for 

L o < x < L o 



y m * ( x _L o - 2R ) - 

for 

L 0 + 2R < x 






m even, antisymmetric modes: 



For this case, the function y m * (z) is defined 

by Equation 2.27 and 



■ y m * ( L o " x ) ■ 

for 

0 < x < L 0 

y m (*) 

= < 

x - L 

r y«*(o> . 

for 

L o < x < L o 



y m * ( x -l 0 -2R) , 

for 

L 0 + 2R < x 


(3.1) 


(3.2) 
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It will be recalled that y m (x) as constructed is defined only up to the multiplicative constant 
C lm . This constant can be uniquely defined for each mode by requiring that the orthonormality 
condition 



(x) dx 


Mp 


(3.3) 


be satisfied for all m. 

Besides predicting the deflection, we can predict the bending moment at each point of the sys- 
tem if we know the curvature of each mode. These are obtainable from the equation 

y*"(z) = A. 2 (C, coshA. z + C 0 sinh A. z - C- cos K z - C 4 sin k z) , (3.4) 

J m x * m \ lm m 2m m 3m m 4m m f ' 

where the constants of integration are identical with those derived for y m (x). We can obtain the curva- 
ture of the mode y m (x) defined over the full system length by applying the construction rules that gen- 
erated y m (x) from y * (z); that is, Equations 3.1 and 3.2. 

4. TRANSIENT AND STEADY-STATE SOLUTIONS 

With a complete set of modes, it is possible to derive a series solution for Equation 1.1— the 
equation of damped forced vibration of the entire system. 

Substituting Equation 1.2 in Equation 1.1 gives 

CO 


n - 1 

Using Equation 1.6 further reduces this to the form 



_d _ 2 
dx ^ 


d 2 y n (x) 


d 2 a (t) 


da n (*) 


a n( t ) ^I EI ( X ) ■' dx 2 + Kx)y n (x) ^2 + /5(x)y n (x) dt 


dt " 


- F(X, t) . 


(4.1) 


00 

L 


d 2 a n (t) da„(t) 

/°(X) y n (X) ^ + 2 plco n y n (x) ^ + p(x ) y n (x)w n 2 a n (t) 


= F(x, t) , 


(4.2) 


where the additional substitution 


/3(x) - 2p& n 


(4.3) 


has also been made. The term £ is commonly referred to as the "damping ratio" and is related to 
the constants of the system by Equation 4.3. It should be noted that the parameter p is the mass 
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per unit length of the beam alone. The equation states that /3(x) is independent of x; hence, energy 
is uniformly dissipated along the entire system length. This is an approximation, since actually 
zero energy is dissipated across the length of the rigid satellite body. 


The summation and all dependence upon the coordinate x can be removed from Equation 4.2, 
by applying the orthonormality condition stated in Equation 1.8. Multiplying Equation 4.2 through 
by y n (x) and integrating over the entire system length reduces this equation to 


d2 3 n CO 

dt 2 


da n ( t ) 

+ 2 £ n "n 3t + “n 2a n (<0 



(x)F(x, t ) dx 


NLj. 


(4.4) 


where £ n is the damping ratio associated with the nth mode and is given by 


l 


n 


i 

m t 



(x) dx . 


(4.5) 


n odd; symmetric modes excited: 


r = J_ f 

- M t 1 


py n 2 (x) dx 


1 - 


2W t 


gM T y n gM T y n (L 0 ) 


(4.6) 


n even; antisymmetric modes excited: 


r 


m t i py n 2 (x)dx 


2W T W s 

1 " y n 2 (°) ■ 3gM T y. 


n 2 (L 0 )] 


(4.7) 


Equation 4.4 is a differential equation that can be solved by means of the Laplace transform. 
Hence, 


a n 0> 



(X)F(X, 


s ) dx 


s 2 + 21 o} s + co 2 

n n 


(4.8) 
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where 


£{F(x, t)} = F(x, s) 

C^n^)} = 3 n ( S ) ^ 4 * 9 ^ 

and the initial values of a n (t ) are assumed equal to zero. 

The system being linear means that we can use the principle of superposition to determine 
the actual system response to the force distribution F(x, t), by adding the response due to each 
individual force distribution of an equivalent force system. 

Let the equivalent force system be defined by the two force distributions F x (x, t) and F 2 (x, t), 
where 

F x (x, t) = F x (x) f(t) = a force distribution of periodic intensity that is symmetric about the 

system's axis of symmetry and excites only the symmetric modes of 
vibration, 

F 2 (x, t) = F 2 (x)f(t) = a force distribution of periodic intensity that is antisymmetric about 

the system's axis of symmetry and excites only the antisymmetric 
modes of vibration. 

The function f(t) appearing above is a continuous function of time; it describes the intensity 
of the force distribution F x (x, t) and F 2 (x, t) at time t. It will be assumed for this analysis that 
f(t) is periodic with period t 2 and expressible as a Fourier series. That is, 


f(t) 


2 b o 



t + d sin 

m 



(4.10) 


where b and d are the Fourier coefficients associated with the function f(t) given by: 

m tn 


b 


m 


d 


m 



2nr77 

f ( t ) cos -j — t dt 


2lT)77 

f ( t ) sin -Z — t dt 
2 


m = 0, 1, 2 • • • 


m = 1, 2, ■ • • 


(4.11) 


(4.12) 


F(x, t) as it appears in Equation 4.8 can be replaced by the equivalent- force system defined 
above. Then, a direct substitution of Equation 4.10 in Equation 4.8 yields an expression that can 
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be readily inverted by convolution: 


*„(»> = 


M T-n (l-«n 2 ) 


1/2 


1 . V -1 / 2m 7T 

2 b o + /_ J (b B c ° s ^7 

m= 1 ' 


2m 77 

+ d sin — z — 

m t j 


•)1 


* e n sin a> n (l - £ n 2 ) 1/2 t 


(4.13) 


where 

n odd; symmetric modes: 


K 


n 



y n ( x ) F i ( x ) dx 


(4.14) 


n even; antisymmetric modes: 


K n “ [ y n ( x ) F 2 ( x ) dx (4.14a) 


and is the symbol used to denote convolution. 

Once evaluated, the expression for a n (t) can be substituted in Equation 1.2 and the displace- 
ment at any point x along the system calculated for any time t. Furthermore, differentiating 
Equation 1.2 twice and multiplying the results by El gives the bending moment at any point for any 
time t. 

Letting 


i 1 ~ £ n 2 ) 1/2 = damped natural frequency of nth mode, 

- £ n w n = inverse of time constant for nth mode, 

0 - 2rr/ 1 2 = frequency of applied- force distribution, 
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and performing the convolution operation called for in Equation 4.13, leads to 


a n( t ) “ M,. 


r ft 

n 11 I -v t 

fi n \2 b ° J o 6 " Sinfi n 


r dr 


5 f- f 


cosm#(t“T)e n sinn n rdr+d m I sinm0(t-T)e n sinO n dr 


r . 

I sir 

Jo 


(4.16) 


K - \ 


n 

n 

e n [v sirtH t + 0 cos 0 t) 

\, n n n n / 

2m t n n 

f* 

v 2 + n 2 

n n 

V 2 +fi 2 J 

n n 


CO 

m - 1 


fa. 


n£>) cos m& t + v n sin mO t e ^ ' [ v n sin n n t + (fl„ ' md) cos t] 


(n n - m ^ 


V 2 + (^n ~ m ^) 2 


(n n + m#) cos m<9 t - v n sin m6 \ 


-v t - 

e n v sin H t + 

n n 


z> n 2 + (Q n + m#) 2 


(0 R + mO) cos tj 
v n 2 + ( 0^ + m<9) 2 


L 


+ 7 d_ 

m = 1 


(n n -md) sinmtf t - v n cos m0 t e ' [ v n cos n „ 1 " ( Q n ~ m5 ) sir,n n l ] 

v n 2 + (n n - me) 2 + v* + (n n - m e) 2 

(n n + me) sin me t + v n cos me t e " ' [ v n cos 1 ' ( n n + ml9 ) sin fi n *] 


v n 2 + (0 n + m#) 2 


V n 2+ ( n n +m ^) 2 


(4.17) 


This equation for a n (t) can be used to define the transient response of the system when ini- 
tially at rest and then set in vibration by the defined applied force F(x, t). In order to investigate 
the response at the point x of the system at time t, we substitute Equation 4.17 in 


Y(x, t) 


00 



n - 1 


to yield the transverse displacement, and in 


El Y" (x, t) 


ei J2 'n(‘)y„'W 

n - 1 


(4.18) 


(4.19) 


to yield the bending moment. 
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By observing how each term 


a n (t)y n (x) (4.20) 

and 

E 1 a n ( t )y n " (x) (4.21) 


changes as a function of time, the contribution to the total response from each individual mode can 
be ascertained for the particular forcing function in question. 

In addition to the transient response, the steady- state response of each mode to the particular 
forcing function is of interest. In deriving the resultant steady- state response of the system, note 
that the steady- state generalized coordinate a n (t) ST given by 


a n ( t )sT 


K 


2M t \ 



+ n 2 


+ 


00 

L 


(0 n - mi^) cos m# t + s in m# t (Q n + m(9) cos mO t - v ^ s in rc\0 t 

V n + ( n n + md Y + y n 2+ ( n n + me ) 2 


00 

L 


(0^ - mtf) s in t - i/ n cos m# t [p- n + m(9) s in mO t + v ^ cos m 0 t 
v 2 + (Q n - m ^) 2 + v 2 + (Q n + m u) 2 


(4.22) 


is not necessarily in phase at the steady- state time t for all n, since it is obtained by evaluating 
a n (t)with the exponential terms deleted. 


5. SOLUTION TO FREQUENCY EQUATION 

It has been shown that the frequency equation associated with the symmetric modes of vibra- 
tion of the entire system is 


W T 


w s 

tanh XL Q + tan \L 0 + 2\v™ (* + sec sec ^ ^o) 

W s 

2XL o ■ ( XL o) 2 ( tan XL o “ tanh ' VL o) 


(5.1) 
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and the frequency equation associated with the antisymmetric modes of vibration of the entire 
system is 


W T 


ft anh \L 0 - tan X.L Q - 2^L 0 R* tan \L 0 tanh A.L 0 - (^L 0 R*) 2 (tanh \L 0 + tan X.L 0 ) 



- 3J (\L 0 ) 3 (l + sec KL q sechXL 0 ) 
2 tan XL 0 tanh X.L Q + 2R* kh Q (tanh + tan \L 0 ) + 2 (R* A.L q ) 2 

+ 3j (\L 0 ) 3 (tanh A.L 0 - tan \L Q 


(5.2) 


By the proper definition of the parameters contained in these equations it is possible to pre~ 
diet the natural frequencies of any system that may be modeled as the symmetric double-beam 
system shown in Figure 1. If the system to be modeled as a double-beam system has four or more 
elastic members attached to a rigid central body, the frequency equations given above yield only 
those frequencies associated with the system mode shapes that are similar to the symmetric and 
antisymmetric modes defined. 


The physical parameters needed to solve the frequency equations of the double beam model 

are: 

W s = weight of rigid satellite body (without beams and tip weights), 

I = mass moment of inertia of rigid satellite body (without beams and tip weights) about the 
axis through the center of mass and normal to the plane of vibration, 

I B = mass moment of inertia of a beam, of length 2L 0 , about an axis normal to the length and 
passing through its centroid, 

W T = tip weight, 

L 0 = length of a single beam, 
p = mass per unit length of beam, 

El = bending stiffness of beam, 

R = radius of symmetric satellite body, 
g = acceleration of gravity. 

Equations 5.1 and 5.2 show that the solutions to the frequency equations will depend only on 
the magnitudes of four distinct dimensionless ratios. 


In particular, the solutions for the symmetric modes involves parameters 


\ ’ W B 


while the solutions for the antisymmetric modes involves 


W B ’ 



R* . 
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These can be calculated for the symmetric double- beam system from the above definitions and 
the equations 

W B = PgL 0 , 


= T^o 3 •, < 5 -3) 

Figures 3 through 7 attempt to provide sufficient information so that Equations 5.1 and 5.2 may be 
accurately solved for all practical values of the defined ratios. With these approximate solutions 
and the relation 



(5.4) 


it is possible to obtain a very good estimate of the first few natural frequencies of the satellite 
system in question. 


i- 

£ 

o 

t— 
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Figure 3— Symmetric mode, W T /W B vs \L 0 for 
different values of W s /W B . 
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Each curve labeled in Figure 3 plots w T/ /w B vs KL 0 as obtained from Equation 5.1 for the value,* 
of W s/ /w B tabulated on the figure. This figure contains sufficient information to give an accurate 
estimate of the first two values of A.L 0 that satisfy Equation 5.1. The natural frequencies associatec 
with the estimated values of \h Q correspond to the first two modes of vibration for which the 
satellite body can translate but not rotate. For the double-boom system, these are the first two 
symmetric modes of vibration. 

A set of curves which provide sufficient information to accurately estimate the natural fre- 
quencies associated with the antisymmetric modes of vibration can be obtained by plotting 
W t/ % vs XL 0 (as derived via Equation 5.2) for incremental values of both l S/ /l B and R* in the 
range 


0 



< 00 


and 


0 < R* < o. 10 . 

Since this involves incrementally varying three parameters; the following approach is used to pre- 
sent the results in a more compact form. 

Intuitively one fells that the least critical parameter appearing in Equation 5.2 is R* . That is, 
for satellites having long elastic members such that 

0 < R* < 0. 1 , 

the natural frequencies will be essentially the same, no matter what the exact value of R* may be. 

Figures 4 and 5 lend graphic support to this intuitive notion and thus provide one with suf- 
ficient information to determine the magnitude range of inertia ratio l s /l B over which R* plays an 
insignificant role in the estimation of the system T s first two antisymmetric modal frequencies. 

That is, we show how sensitive the asymptotes and zero crossings of Equation 5.2 are to change in 
i s / i b and R * • Since it will be shown that the curves W T /W B vs a.l 0 as derived from Equation 5.2 
have shapes similar to those shown in Figure 3, it may be inferred that any intermediate value of 
W T /W B between zero and infinity will exhibit the same sensitivity to changes in l s /l B and R* as do 
the asymptotes and zero crossings. 

Figure 4 plots I s /l B vs XL 0 for the values of R* tabulated on the figure. These curves are ob- 
tained from Equation 5.2 when W T /W B is set equal to infinity; that is, from 

j 6 [^tan\L 0 tanh\L 0 +R*\L 0 (t anh \L Q + tan kh Q ) + (R* AL 0 ) 2 ] 

(^L 0 ) 3 (tan \L 0 - tanh XL 0 ) (5*5) 
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SATELLITE -TO -BEAM -INERTIA L/I. SATELLITE -TO -BEAM- 1 NERTA RATIO l s /l 

^ r- 1 



XL Q ( radians) 


Figure 4— Antisymmetric mode, l s /l„ vs\L 0 for different 
values of R*, when W T /vV B = 



Figure 5— Antisymmetric mode, l s /l B vs XL Q for different 
values of R*, when W T /W B “ 0. 



■ i am linn uni i hi 11 


■ »•■ 


1 1 


Similarly Figure 5 plots l s /l B vs KL 0 for the same values of R*, as obtained from Equation 5.2 when 
W T/ /W B is set equal to zero; that is, from 


I 

I 


s 

B 


|\anh XL 0 - tan A.L Q 


(5.6) 


Making use of the fact that the solutions to Equation 5.2 are virtually insensitive to changes in 
the parameter R* over the limits defined, Figures 6 and 7 plot W T /w B vs A.L 0 as derived from Equa- 
tion 5.2 for the values of R* = 0.01 and R* = 0.10 respectively and the values of I s /l B listed on 
the figure. These figures contain sufficient information to give accurate estimates of the first two 
values of XL 0 that satisfy Equation 5.2 for any system that can be modeled as discussed. The 
natural frequencies associated with the estimated values of XL 0 correspond to the first two modes 
of vibration for which the satellite body can rotate but not translate. For the double-beam system 
these are the first two antisymmetric modes of vibration. 

Aside from providing the information needed for solving both frequency equations, Figures 3 
through 7 along with Equation 5.4 provide quantitative and qualitative insight as to how changes in 
the system parameters will change the system’s natural frequencies. 

The range of W s /w B values listed on Figure 3 shows that the root boundary conditions for the 
symmetric modes vary between zero satellite weight (guided root) and infinite satellite weight 
(clamped root). Similarly, the range of the ratio magnitudes listed on Figures 6 and 7 show that 
the root boundary conditions for the antisymmetric modes vary between zero satellite rotational 
inertia (hinged root) to infinite rotational inertia (clamped root) for each of the listed ratios R*. 

For both cases, the variation in the tip boundary conditions from zero tip weight (free tip) to in- 
finite tip weight (hinged tip) is illustrated by choosing W T /W B as the vertical axis of Figures 3, 6, 
and 7. 

Figure 8 gives a ready reference for the general solution of the frequency equations in the 
limiting cases mentioned above. The figure shows the schematic diagram of each limiting case 
along with the solutions to its frequency equation. Reference 3 gives the solutions for these cases, 
also for cases corresponding to other common beam boundary conditions. 

Figures 3 through 7 and Equation 5.4 suggest the following observations: 

1. The solutions VL 0 of the frequency equations are independent of the bending stiffness El . 
Thus the system’s natural frequencies are directly proportional to the square root of the bending 
stiffness. 

2. The solutions A.L 0 of the frequency equations have nonlinear dependence on the boom weight 
W B . Thus a simple proportionality relation between boom length L 0 or mass distribution p and the 
natural frequencies exists only in special cases. 
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TIP-TO-BEAM-WEIGHT RATIO W T /W 
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Figure 7— Antisymmetric 
W T /W B vs \L 0 for d 
values of l s /l B , when R’ 
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Figure 8— Natural frequencies of uniform beams with 
different boundary conditions. 


3. As the satellite-to-beam-weight ratio 
W s/ /w B increases, the solution (^L 0 ) x associated 
with the first symmetric mode of vibration ap- 
proaches the value that results from a clamped- 
root condition. Figure 3 shows that for W s/ /w B 
greater than 100 and W T/ /w B less than 20 the root 
may be assumed clamped and the frequency 
associated with the first symmetric mode of 
vibration is approximately that of a clamped 
beam with tip weight. 

4. Figures 4 and 5 show that for satellites 
having a satellite-to-beam-length ratio R* in 
the interval 


0 < R* < 0. 10 , 

the solution (>JL 0 ) 2 associated with the first 
antisymmetric mode of vibration will be es- 
sentially independent of the magnitude of R* 
for all inertia ratios l s /l B greater than 0.10 
For the range 


T 1 < 0. 10 


and 


0.05 < R* < 0. 10 , 


the quantity R* becomes a significant parameter in Equation 5.2. This situation implies that the 
satellited rotational inertia is small compared with the boom’s; at the same time, the boom cannot 
be extremely long compared with the length of the rigid satellite body (a situation approaching that 
of a physically unrealistic system). Thus, for most actual satellites, the magnitude of R* and 
hence the actual satellite body size will be of second order importance in predictions of the sys- 
tem’s first antisymmetric modal frequency. 

5. Figures 6 and 7 show that as the satellite inertia-to-beam inertia ratio i s /i B increases, 
the solution (al 0 ) 2 associated with the first antisymmetric mode of vibration approaches the value 
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that results from the clamped-root condition defined by curve number 11 on both figures. The in- 
ertia ratio needed to effectively clamp the root and hence restrict the rotational motion of the 
satellite body is extremely large and strongly dependent on the tip-to-beam-weight ratio w T/ /w B . 
For example, if W T/ /w B = 1, then the inertia ratio l s /l B must be greater than 50 for the root to be 
considered clamped. 


6. Comparing the values of (A.L 0 ) 2 that define the asymptotes of the curves shown in Figures 
6 and 7 with the true values of the asymptotes given in Figure 4 shows that for every satellite 
having an inertia ratia l s /i B less than 5, a tip-to-beam-weight ratio can be found such that the 
system will vibrate as if the system had infinitely heavy tip weights. For example, if the inertia 
ratio of the satellite were less than 1, then the system would vibrate as if its tips were pinned in 
inertial space by infinitely heavy tip weights, for any W T/ /W B greater than 5. 

7. As w T /w B increases from 0 to 3, the solutions (\L 0 ) n , n = 3, 4, • • • of the frequency equa- 
tions, associated with the second and higher symmetric and antisymmetric modes, undergo a rapid 
decrease irrespective of satellite body weight, size, and rotational inertia. As W T /W B is increased 
beyond 3, (\L 0 ) n becomes insensitive to any further change and approaches the solution associated 
with the case of a hinged tip. Thus, if the tip-to-beam-weight ratio exceeds 3, the symmetric 
double-beam system will have its tips effectively hinged in inert^JL space for the third and all 
higher modes of vibration. 


6. DYNAMIC RESPONSE OF ATS-D 

The introduction stated that this paper would solve a particular problem: to determine 
whether the microthruster on ATS-D could excite a resonant response of sufficient magnitude to 
buckle the long elastic members (booms) attached to the spacecraft body. 

In order to introduce conservatism into the model and reduce the problem to a form com- 
patible with this analysis, we make some simplifying assumptions: 


1. The actual satellite deploying four booms with tip weights may be adequately approximated 
by the symmetric double-beam system. 

2. The microthruster acts normal to the system at the point on the satellite body coincident 
with the point at which one of the beams is clamped. The intensity of the microthruster may be 
approximated by a periodic triangular wave having a rise time unequal to the decay time and a 
peak intensity of one micropound. 
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Let 

t x = rise time of triangular wave, sec, 
t = period of triangular wave, sec, 

10' 6 = one micropound, lb, 

8(x) = Dirac Delta function. 

Then the actual applied force distribution F(x, t) is described by 

F(x, t) = $(x-L- R) f(t) , (6.1) 


where, in the time interval 


0 < t < t 2 , 


f(t) = 10' 6 for 0 < t < tj 


f(t) = 10 


H’- tj f ° r 


<*^2- 


( 6 . 2 ) 


This force distribution is replaced in the analysis by the equivalent force system described by 
Fj (x. t) and F 2 (x, t); that is: 


Fj(x,t) = Fj (x) f ( t ) = S(L-X)f(t) 

F 2 (x,t) = F 2 (x)f(t) 


( R \ / 

R \" 

S( L +~ 2 ~ xl - Sll 

_ T“ xj 


f(t). 


(6.3) 

(6.4) 


where 


F t (x, t) = force applied at the center of mass, 

F 2 (x, t) = couple applied about the center of mass. 

It follows from Equations 4.11, 4.12, and 6.2 that the Fourier coefficients described in Section 
4 are: 

b 0 - 10-* 


b 


d 


m 


10 " 6 1 2 2 

2m 2 7T2t 1 (t 2 -t ] )_ 
10' 6 t 2 2 

2m 2 v 2 ti(t 2 -tj) 


1 - COS 


2m7 t 


2mn 

XT 


i 



m = 1,2 (6.5) 
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and that the constants K n described by Equations 4.14, 4.14a, 6.3, and 6.4 are: 

n odd; symmetric modes: 

K n = y n (L) , (6.6) 

n even; antisymmetric modes: 

K n = y n ( L + f) - (l - f) • ( 6 - 7 > 

3. The vibrational motion will be damped by the internal friction of the boom material only 
and this is describable by a viscous-type damping term in the analysis. Reference 4 presents the 
results of experiments aimed at precisily defining the magnitude of the damping coefficient £ for 
De Havilland-type booms. It may be concluded from this that worst- case- type results may be ob- 
tained by setting 

£ = 1CT 4 . (6.8) 

4. The adjusted geometric and physical parameters that describe the double -beam model of the 
four boom ATS-D satellite are assumed to be approximated by the magnitudes listed in Table 1. It 
must be noted that the mass distribution, bending stiffness, and tip weight of the beam are exactly 
twice the corresponding values for a single boom of the actual satellite system. 


Table 1 

Adjusted Parameters for Double -Beam Model of ATS-D. 


Symbol 

Parameter 

Value 

P 

Mass distribution of boom 

0.8125 x 10 ~ 3 slugs ft 1 

El 

Bending stiffness of boom 

26 lbs ft 2 

L„ 

Boom length 

123.5 ft 

R 

Radius of ATS 

2 ft 

W T 

Tip weight 

16 lb 

w s 

Weight of ATS (without booms and tip weights) 

641.4 lb 

K 

Rotational inertia of ATS (without booms and tip weights) 

76.9 slugs ft 2 


To give accurate values of the first n natural frequencies of the double-beam system, Equa- 
tions 5.1 and 5.2 must be solved numerically. It will be recalled that Equation 5.1 is the frequency 
equation for the symmetric modes and that Equation 5.2 is the frequency equation for the antisym- 
metric modes. 

Table 2 lists the first ten values of \L 0 that satisfy these equations. The odd-numbered values 
correspond to the solutions to Equation 5.1, and the even-numbered values correspond to the solutions 
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Table 2 

First Ten Natural Frequencies of ATS-D. 


K)n 

(radian) 

CO 

n 

(radian sec -1 ) 


T 

n 

(sec) 

K), 

0.88243 


0.009133 

T i 

687.99 

o 

lO 
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Figure 9 — First symmetric mode of ATS system 
described in Table 1 . 


to Equation 5.2. The table also lists the modal 
frequency and period associated with each solu- 
tion (^L 0 ) n . 

Both the symmetric and antisymmetric 
mode shapes of the entire system can be gen- 
erated by using the values of (^L 0 ) n listed in 
Table 2 and the other quantities given in Tablet. 

Figures 9 through 13 plot the normalized 
displacement y n (x) vs the position coordinate x 
for the first five symmetric modes. The equa- 
tions used to define these mode shapes are ob- 
tained by evaluating Equation 2. 23 for those Ml 0 T s 
with odd subscripts listed in Table 2, applying 
the construction rules given by Equations 3.1, 
and normalizing according to Equation 3.3. 



Figure 10— Second symmetric mode of ATS system 
described in Table 1. 



described in Table 1. described in Table 1. 
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Figures 14 through 18 plot the normalized 
displacement y n (x) vs position coordinate x for 
the first five antisymmetric modes. The equa- 
tions used to define these mode shapes are ob- 
tained by evaluating Equation 2.27 for those 
\L 0 T s with even subscripts listed in Table 2, ap- 
plying the construction rules given by Equation 
3.2, and normalizing according to Equation 3.3. 

Besides the mode shapes of the system, it 
is also possible to generate equations that de- 
fine the normalized curvature of every point 
along each mode shape. This may be done by 



Figure 14— First antisymmetric mode of ATS system 
described in Table 1. 



Figure 13— Fifth symmetric mode of ATS system 
described in Table 1. 



Figure 15 — Second antisymmetric mode of ATS system 
described in Table 1. 



Figure 16— Third antisymmetric mode of ATS system 
described in Table 1 . 


Figure 17— Fourth antisymmetric mode of ATS system 
described in Table 1. 
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substituting the values of c ln determined from the normalization of y n (x) in Equation 2.21 or 2.25 
to determine C 2n , C 3n , and C 4n . These normalized constants can then be substituted in Equation 
3.4 to give an expression for y n " (z). This resulting expression can be substituted in Equation 3.1 



Figure 1 8 — Fifth antisymmetric mode of ATS system 
described in Table 1. 


or 3.2, depending on whether n is odd or even, 
respectively, to construct y n " (x), the normalized 
curvature of the nth mode. 

Figures 19 through 23 plot the normalized 
curvature y n " (x) of the first five symmetric 
modes of the system (n odd) vs position coor- 
dinate x. 

Figures 24 through 28 plot the normalized 
curvature y n " (x) of the first five antisymmetric 
modes of the system (n even) vs position coor- 
dinate X. 



Figure 20— Curvature of second symmetric mode. 

Figure 19— Curvature of first symmetric mode. 




Figure 21— Curvature of third symmetric mode. 


Figure 22— Curvature of fourth symmetric mode. 






From Figures 9 through 28, the following observations can be made concerning the possible 
effects of vibrational motion on the ATS-D system: 

1. Figures 9 through 18 show that the satellite body and tip weights will be essentially pinned 
in inertial space for all modes of vibration except the first symmetric mode, 

2. Figures 14, 15, and 16 show that if any one of the first three antisymmetric modes of vibra- 
tion are excited, the rotational inertia of the satellite body is not large enough to prevent significant 
rotational motion of the body about the center of mass. 

3. For all modes of oscillation except the first antisymmetric mode, the peak curvature and 
hence the peak bending moment will exist at the boom's root. For the first antisymmetric mode 
(see Figure 24) the point of peak curvature will occur at a point 45 feet from the root. The curva- 
ture at this point will be about 1.5 times the curvature at the root. 

The object of this example is to determine whether the viscoelastic damping of the boom ma- 
terial itself can prevent the bending moment at any point along the boom from exceeding its criti- 
cal buckling moment when the system is excited at a resonant frequency. 

This may be done by directly comparing the critical buckling moment and peak bending moment 
of the system for each excitation frequency under study. We can avoid evaluating the complete ex- 
pression defining the peak bending moment, 


max max El Y" (x, t) - max max El / a ( t ) y " (x) , (6.9) 

x t x t / > n n 

n = 1 

by using the following fact: when the system is excited by the triangular wave whose period equals 
one of the natural periods of the system, the major contribution to the resultant bending moment at 
any point is associated with the particular mode excited at resonance. 

Thus, when -the excitation period t 2 defined in Equation 6.2 is given by 


and the rise time t x is given by 


2 77 
CO 

n 


3 


2 n 

3-n 


the peak bending moment is approximated by the expression 

max max El Y" (x, t) = max max El a ( t ) y n " (x) . 

X t x t 


( 6 . 10 ) 


( 6 . 11 ) 


( 6 . 12 ) 
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It will be recalled that for all modes of vibration except the first antisymmetric (n = 2) the 
peak curvature is at the boom root; i.e., at x = L±R, and that 

|y n *(L+R)| = ly n "(L-R)| . (6.13) 

Therefore, by applying this fact the maximization on x called for in Equation 6.12 can be per- 
formed. This leads to the expression 


max max El Y" (x, t) 

X t 


pi max a n ( t ) y n " (L + R) , 

for 

n / 2 

[l. 5 El max a n (t)y n " (L + R) , 

for 

n = 2 


(6.14) 


where, for n = 2, the peak curvature is 1.5 times 
the curvature at the root. The maximization on 
t called for in Equation 6.14 can be performed 
by directly substituting Equations 6.3 through 
6.8 in Equation 4.22 and computing the maxi- 
mum value of a (t) CT in the interval 0 < t < t_. 

The points marked on Figure 29 show what 
the contribution to the peak bending moment at 
the boom root would be for each mode if it were 
excited at its natural period. That is, for the 
excitation frequency 


0.1 

n 


277 


(6.15) 


shown on the horizontal scale, the magnitude of 


max max El a n ( t ) y n " ( x ) 


can be ascertained from the vertical scale. 
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Figure 29— Peak bending moment contribution at boom 
root from each mode for a triangular forcing function 
having a frequency equal to natural frequency. 


Thus, even if the damping ratio were as low as 10 4 and the microthruster did excite a reso- 
nant frequency, the critical buckling moment— which is about one ft -lb— would not be exceeded. 


SUMMARY 

The equations developed herein may be applied to any satellite having long elastic members 
that can be modeled as a symmetric double-beam system, as shown in Figure 1. The following is 
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a summary of the equations and the conclusions that can be drawn by a close examination of them 
and the figures. 

1. Because of the system’s symmetry, the normalized modes of vibration can be separated 
into two groups: symmetric and antisymmetric, respectively, about the system’s axis of symmetry. 
Equations 5.1 and 5.2 are the frequency equations of the symmetric and antisymmetric modes, 
respectively. 

2. The frequency equations are written in terms of the dimensionless parameters: 

(a) w T /W B = ratio of tip to beam weight 

(b) r/l 0 = ratio of satellite body radius to beam length 

(c) W s /\V B = ratio of satellite body to beam weight 

(d) I S /I B = ratio of satellite body to beam inertia 

The natural frequencies of the system may be determined with slide-rule accuracy by simply 
defining the above parameters and interpolating between the appropriate curves on Figures 3, 6 
and 7. 

3. The following conclusions may be drawn from Equations 5.1 and 5,2 and the corresponding 
figures. 

(a) For W s /W B >100 and W T/ /W B < 20, the frequency associated with the first symmetric mode 
of vibration is approximately that of a clamped beam with tip weight. 

(b) For I s/ /l B >0.1, 0 <R/L 0 < 0.1 and 0 <W T/ /W B < <», the frequency associated with the first 
antisymmetric mode of vibration is approximately independent of R/L 0 . 

(c) For I s /l B < 0.1, 0.05 <R/L 0 <0.1, and 0 £W T /W B <<», the quantity R/L 0 becomes a signifi- 
cant parameter in the prediction of the first antisymmetric modal frequency. 

(d) For l s/ /i B > 250 and W T/ /W B < 30, the frequency associated with the first antisymmetric 
mode of vibration is approximately that of a clamped beam with tip weight. 

(e) For I s /I B < 5 and W T/ /w B >10, the frequency associated with the first antisymmetric 
mode of vibration is approximately that of a similar system having its tips hinged in 
inertial space. 

(f) For w T/ /w B > 3, all higher symmetric and antisymmetric modes of vibration can be ap- 
proximated by a similar system having its tips hinged in inertial space. 

4. The equations defining the symmetric and antisymmetric mode shapes are developed in 
terms of dimensionless parameters and are given by Equations 3.1 and 3.2. 
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5. The equations defining both the transient and the steady- state response of the system in 
terms of its normal modes of vibration are derived for an arbitrary force distribution of periodic 
intensity. These Equations are 4.17, 4.18, and 4.22. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, February 5, 1968 
630-12-02-01-51 
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Appendix A 


Symbol List 

a n (t) time- dependent generalized displacement coordinate associated with the nth mode, ft. 

b m mth Fourier cosine coefficient. 

c i* C 2 , c 3 , C 4 constants of integration, 
d mth Fourier sine coefficient. 

m 

El(x) bending stiffness at x, lb ft 2 . 

El bending stiffness along length of elastic beam, lb ft 2 . 

F( x, t) force distribution at (x, t), lb ft -1 . 

F i ( x > t ) F i ( x ) i ( t ) force distribution of periodic intensity that is symmetric about the sys- 
tem^ axis of symmetry and excites only the symmetric modes of vibration, lb ft -1 . 

F 2 (x, t) f 2 (x) f(t) force distribution of periodic intensity that is antisymmetric about the 
system T s axis of symmetry and excites only the antisymmetric modes of vibration, 
lb ft’ 1 . 

g acceleration of gravity, ft sec -2 . 

I B mass moment of inertia of a beam of length 2L 0 about an axis normal to the length 

arid passing through its centroid, slugs ft 2 . 

I s mass moment of inertia of rigid satellite body (without beams and tip weights) about 

the axis through the center of mass and normal to the plane of vibration, slugs ft 2 . 

2L length of entire system, ft. 

L 0 length of one beam, ft. 

M t total mass of entire system, slugs. 

m, n integers 0, 1, 2 • • • 

2R total length of rigid satellite body, ft. 

2R * normalized length of rigid satellite body. 
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t time, sec. 

tj rise time of triangular wave, sec. 

t 2 period of triangular wave, sec. 

W B weight of single beam without tip weight. 

W T tip weight, lb. 

W s weight of rigid satellite body (without beams and tip weights), lb. 

w; effective inertia weight of satellite body, lb. 

x position coordinate measured positive from the beam tip, ft. 

y n (x) nth normalized mode of the undamped unforced equations of motion of the entire 

system. 

y n * (z) nth mode of the undamped unforced equation of motion of one beam, used to con- 

struct y n (x). 

Y(x, t) transverse displacement at (x, t), ft. 

z position coordinate defined over flexible portion of one beam, ft. 

/3(x) viscous damping coefficient at x, lb sec ft~ 2 . 

J3 viscous damping coefficient of the beams alone, lb sec ft 2 . 

S(x) Dirac Delta function. 

Sn, m Kronecker Delta function. 

P(x) mass distribution of the satellite beam system at x, slugs ft" 1 . 

p mass distribution along the beam, slugs ft" 1 . 

E n damping ratio associated with the nth mode. 

£ damping ratio measured experimentally. 

nth eigenvalue of the defined eigenvalue problem. 

<x> n nth natural frequency of system. 

0 frequency of applied force distribution. 

v n inverse of the time constant associated with the nth mode of vibration. 

H n damped natural frequency of the th mode of vibration. 
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